CRMnet: A deep learning model for predicting gene expression from large regulatory sequence datasets.

Recent large datasets measuring the gene expression of millions of possible gene promoter sequences provide a resource to design and train optimized deep neural network architectures to predict expression from sequences. High predictive performance due to the modeling of dependencies within and between regulatory sequences is an enabler for biological discoveries in gene regulation through model interpretation techniques. To understand the regulatory code that delineates gene expression, we have designed a novel deep-learning model (CRMnet) to predict gene expression in Saccharomyces cerevisiae. Our model outperforms the current benchmark models and achieves a Pearson correlation coefficient of 0.971 and a mean squared error of 3.200. Interpretation of informative genomic regions determined from model saliency maps, and overlapping the saliency maps with known yeast motifs, supports that our model can successfully locate the binding sites of transcription factors that actively modulate gene expression. We compare our model's training times on a large compute cluster with GPUs and Google TPUs to indicate practical training times on similar datasets.

Purification and membrane interactions of human KCNQ1100-370 potassium ion channel.

KCNQ1 (Kv7.1 or KvLQT1) is a voltage-gated potassium ion channel that is involved in the ventricular repolarization following an action potential in the heart. It forms a complex with KCNE1 in the heart and is the pore forming subunit of slow delayed rectifier potassium current (Iks). Mutations in KCNQ1, leading to a dysfunctional channel or loss of activity have been implicated in a cardiac disorder, long QT syndrome. In this study, we report the overexpression, purification, biochemical characterization of human KCNQ1100-370, and lipid bilayer dynamics upon interaction with KCNQ1100-370. The recombinant human KCNQ1 was expressed in Escherichia coli and purified into n-dodecylphosphocholine (DPC) micelles. The purified KCNQ1100-370 was biochemically characterized by SDS-PAGE electrophoresis, western blot and nano-LC-MS/MS to confirm the identity of the protein. Circular dichroism (CD) spectroscopy was utilized to confirm the secondary structure of purified protein in vesicles. Furthermore, 31P and 2H solid-state NMR spectroscopy in DPPC/POPC/POPG vesicles (MLVs) indicated a direct interaction between KCNQ100-370 and the phospholipid head groups. Finally, a visual inspection of KCNQ1100-370 incorporated into MLVs was confirmed by transmission electron microscopy (TEM). The findings of this study provide avenues for future structural studies of the human KCNQ1 ion channel to have an in depth understanding of its structure-function relationship.

Plasma proteomic analysis of association between atrial fibrillation, coronary microvascular disease and heart failure.

The clinical association between atrial fibrillation (AF), coronary microvascular disease (CMD) and heart failure with preserved ejection fraction (HFpEF) is highly prevalent, however the mechanism behind this association is not known. We hypothesized that plasma proteomic analysis can identify novel biomarkers and the mechanistic pathways in concomitant AF, CMD and HFpEF. To discover circulating biomarkers for the association between AF, CMD and HFpEF, an unbiased label-free quantitative proteomics approach was used in plasma derived from patients who underwent coronary physiology studies (n=18). Circulating proteins were analyzed by liquid chromatography-mass spectrometry and screened to determine candidate biomarkers of the concomitant AF, CMD and HFpEF. We identified 130 dysregulated proteins across the groups with the independent patient replicates. Among those, 35 proteins were candidate biomarkers of the association between AF, CMD and HFpEF. We found significantly elevated SAA1, LRG1 and APOC3 proteins in the coexistence of AF, CMD and HFpEF, whereas LCP1, PON1 and C1S were markedly downregulated in their associations. AF was associated with reduced LCP1, KLKB1 and C4A in these patients. Combined downregulation of PON1 and C1S was a marker of concurrent HFpEF and CMD. PON1 was associated with HFpEF while C1S was a marker of CMD. These proteins are related to inflammation, extra cellular remodeling, oxidative stress, and coagulation. In conclusion, plasma proteomic profile provides biomarkers and mechanistic insight into the association of AF, CMD and HFpEF. SAA1, LRG1, APOC3, LCP1, PON1 and C1S are candidate markers for the risk stratification of their associations and potential underlying mechanistic pathways.

Identification of Alternative Splicing and Polyadenylation in RNA-seq Data.

As well as the typical analysis of RNA-Seq to measure differential gene expression (DGE) across experimental/biological conditions, RNA-seq data can also be utilized to explore other complex regulatory mechanisms at the exon level. Alternative splicing and polyadenylation play a crucial role in the functional diversity of a gene by generating different isoforms to regulate gene expression at the post-transcriptional level, and limiting analyses to the whole gene level can miss this important regulatory layer. Here, we demonstrate detailed step by step analyses for identification and visualization of differential exon and polyadenylation site usage across conditions, using Bioconductor and other packages and functions, including DEXSeq, diffSplice from the Limma package, and rMATS.

A unique histone 3 lysine 14 chromatin signature underlies tissue-specific gene regulation.

Organismal development and cell differentiation critically depend on chromatin state transitions. However, certain developmentally regulated genes lack histone 3 lysine 9 and 27 acetylation (H3K9ac and H3K27ac, respectively) and histone 3 lysine 4 (H3K4) methylation, histone modifications common to most active genes. Here we describe a chromatin state featuring unique histone 3 lysine 14 acetylation (H3K14ac) peaks in key tissue-specific genes in Drosophila and human cells. Replacing H3K14 in Drosophila demonstrates that H3K14 is essential for expression of genes devoid of canonical histone modifications in the embryonic gut and larval wing imaginal disc, causing lethality and defective wing patterning. We find that the SWI/SNF protein Brahma (Brm) recognizes H3K14ac, that brm acts in the same genetic pathway as H3K14R, and that chromatin accessibility at H3K14ac-unique genes is decreased in H3K14R mutants. Our results show that acetylation of a single lysine is essential at genes devoid of canonical histone marks and uncover an important requirement for H3K14 in tissue-specific gene regulation.

Activation of AMP-Activated Protein Kinases Prevents Atrial Fibrillation.

Atrial fibrillation (AF) is common, yet there is no preventive therapy for AF. We tested the efficacy of AMP-activated protein kinase (AMPK) activators, metformin, and aspirin, in primary prevention of AF in cardiac-specific liver kinase B1 (LKB1) knockout (KO) mouse model of AF. Incidence of spontaneous AF was significantly reduced in treated KO mice with metformin (10 mg/kg/day) (8.3% in male and 10.3% in female) and aspirin (20 mg/kg/day) (29.4% in male and 21.4% in female) compared with untreated littermates (81% in male and 67% in female) at 8 weeks (p < 0.05). Prevention of AF was associated with activation of AMPK in treated mice and thereby improvement of mitochondrial function, gap junction proteins (connexin 40/43), and intra- and inter-cellular ultrastructure in atrial myocardium. Fibrosis was significantly less in treated mice atria. Pharmacological activation of AMPK is an effective upstream therapy for the primary prevention of AF in susceptible heart. Graphical abstract.

Active S2168 and inactive S21IRS pinholin interact differently with the lipid bilayer: A 31P and 2H solid state NMR study.

Pinholins are a family of lytic membrane proteins responsible for the lysis of the cytosolic membrane in host cells of double stranded DNA bacteriophages. Protein-lipid interactions have been shown to influence membrane protein topology as well as its function. This work investigated the interactions of pinholin with the phospholipid bilayer while in active and inactive confirmations to elucidate the different interactions the two forms have with the bilayer. Pinholin incorporated into deuterated DMPC-d54 lipid bilayers, along with 31P and 2H solid state NMR (SS-NMR) spectroscopy were used to probe the protein-lipid interactions with the phosphorus head group at the surface of the bilayer while interactions with the 2H nuclei were used to study the hydrophobic core. A comparison of the 31P chemical shift anisotropy (CSA) values of the active S2168 pinholin and inactive S21IRS pinholin indicated stronger head group interactions for the pinholin in its active form when compared to that of the inactive form supporting the model of a partially externalized peripheral transmembrane domain (TMD) of the active S2168 instead of complete externalized TMD1 as suggested by Ahammad et al. JPC B 2019. The 2H quadrupolar splitting analysis showed a decrease in spectral width for both forms of the pinholin when compared to the empty bilayers at all temperatures. In this case the decrease in the spectral width of the inactive S21IRS form of the pinholin showed stronger interactions with the acyl chains of the bilayer. The presence of the inactive form's additional TMD within the membrane was supported by the loss of peak resolution observed in the 2H NMR spectra.

Characterization of the Human KCNQ1 Voltage Sensing Domain (VSD) in Lipodisq Nanoparticles for Electron Paramagnetic Resonance (EPR) Spectroscopic Studies of Membrane Proteins.

Membrane proteins are responsible for conducting essential biological functions that are necessary for the survival of living organisms. In spite of their physiological importance, limited structural information is currently available as a result of challenges in applying biophysical techniques for studying these protein systems. Electron paramagnetic resonance (EPR) spectroscopy is a very powerful technique to study the structural and dynamic properties of membrane proteins. However, the application of EPR spectroscopy to membrane proteins in a native membrane-bound state is extremely challenging due to the complexity observed in inhomogeneity sample preparation and the dynamic motion of the spin label. Detergent micelles are very popular membrane mimetics for membrane proteins due to their smaller size and homogeneity, providing high-resolution structure analysis by solution NMR spectroscopy. However, it is important to test whether the protein structure in a micelle environment is the same as that of its membrane-bound state. Lipodisq nanoparticles or styrene-maleic acid copolymer-lipid nanoparticles (SMALPs) have been introduced as a potentially good membrane-mimetic system for structural studies of membrane proteins. Recently, we reported on the EPR characterization of the KCNE1 membrane protein having a single transmembrane incorporated into lipodisq nanoparticles. In this work, lipodisq nanoparticles were used as a membrane mimic system for probing the structural and dynamic properties of the more complicated membrane protein system human KCNQ1 voltage sensing domain (Q1-VSD) having four transmembrane helices using site-directed spin-labeling EPR spectroscopy. Characterization of spin-labeled Q1-VSD incorporated into lipodisq nanoparticles was carried out using CW-EPR spectral line shape analysis and pulsed EPR double-electron electron resonance (DEER) measurements. The CW-EPR spectra indicate an increase in spectral line broadening with the addition of the styrene-maleic acid (SMA) polymer which approaches close to the rigid limit providing a homogeneous stabilization of the protein-lipid complex. Similarly, EPR DEER measurements indicated a superior quality of distance measurement with an increase in the phase memory time (Tm) values upon incorporation of the sample into lipodisq nanoparticles when compared to proteoliposomes. These results are consistent with the solution NMR structural studies on the Q1-VSD. This study will be beneficial for researchers working on investigating the structural and dynamic properties of more complicated membrane protein systems using lipodisq nanoparticles.

The membrane protein KCNQ1 potassium ion channel: Functional diversity and current structural insights.

BACKGROUND: Ion channels play crucial roles in cellular biology, physiology, and communication including sensory perception. Voltage-gated potassium (Kv) channels execute their function by sensor activation, pore-coupling, and pore opening leading to K+ conductance. SCOPE OF REVIEW: This review focuses on a voltage-gated K+ ion channel KCNQ1 (Kv 7.1). Firstly, discussing its positioning in the human ion chanome, and the role of KCNQ1 in the multitude of cellular processes. Next, we discuss the overall channel architecture and current structural insights on KCNQ1. Finally, the gating mechanism involving members of the KCNE family and its interaction with non-KCNE partners. MAJOR CONCLUSIONS: KCNQ1 executes its important physiological functions via interacting with KCNE1 and non-KCNE1 proteins/molecules: calmodulin, PIP2, PKA. Although, KCNQ1 has been studied in great detail, several aspects of the channel structure and function still remain unexplored. This review emphasizes the structural and biophysical studies of KCNQ1, its interaction with KCNE1 and non-KCNE1 proteins and focuses on several seminal findings showing the role of VSD and the pore domain in the channel activation and gating properties. GENERAL SIGNIFICANCE: KCNQ1 mutations can result in channel defects and lead to several diseases including atrial fibrillation and long QT syndrome. Therefore, a thorough structure-function understanding of this channel complex is essential to understand its role in both normal and disease biology. Moreover, unraveling the molecular mechanisms underlying the regulation of this channel complex will help to find therapeutic strategies for several diseases.

Probing the Dynamics and Structural Topology of the Reconstituted Human KCNQ1 Voltage Sensor Domain (Q1-VSD) in Lipid Bilayers Using Electron Paramagnetic Resonance Spectroscopy.

KCNQ1 (Kv7.1 or KvLQT1) is a potassium ion channel protein found in the heart, ear, and other tissues. In complex with the KCNE1 accessory protein, it plays a role during the repolarization phase of the cardiac action potential. Mutations in the channel have been associated with several diseases, including congenital deafness and long QT syndrome. Nuclear magnetic resonance (NMR) structural studies in detergent micelles and a cryo-electron microscopy structure of KCNQ1 from Xenopus laevis have shown that the voltage sensor domain (Q1-VSD) of the channel has four transmembrane helices, S1-S4, being overall structurally similar with other VSDs. In this study, we describe a reliable method for the reconstitution of Q1-VSD into (POPC/POPG) lipid bilayer vesicles. Site-directed spin labeling electron paramagnetic resonance spectroscopy was used to probe the structural dynamics and topology of several residues of Q1-VSD in POPC/POPG lipid bilayer vesicles. Several mutants were probed to determine their location and corresponding immersion depth (in angstroms) with respect to the membrane. The dynamics of the bilayer vesicles upon incorporation of Q1-VSD were studied using 31P solid-state NMR spectroscopy by varying the protein:lipid molar ratios confirming the interaction of the protein with the bilayer vesicles. Circular dichroism spectroscopic data showed that the α-helical content of Q1-VSD is higher for the protein reconstituted in vesicles than in previous studies using DPC detergent micelles. This study provides insight into the structural topology and dynamics of Q1-VSD reconstituted in a lipid bilayer environment, forming the basis for more advanced structural and functional studies.

Genome-based taxonomic classification.

Bacterial populations are routinely characterized based on microscopic examination, colony formation, and biochemical tests. However, in the recent past, bacterial identification, classification, and nomenclature have been strongly influenced by genome sequence information. Advances in bioinformatics and growth in genome databases has placed genome-based metadata analysis in the hands of researchers who will require taxonomic experience to resolve intricacies. To achieve this, different tools are now available to quantitatively measure genome relatedness within members of the same species, and genome-wide average nucleotide identity (gANI) is one such reliable tool to measure genome similarity. A genome assembly with a gANI score of <95% at the intraspecies level is generally considered indicative of a separate species. In this study, we have analysed 300 whole-genome sequences belonging to 26 different bacterial species available in the NCBI Genome database and calculated their similarity at the intraspecies level based on gANI score. At the intraspecies level, nine bacterial species showed less than 90% gANI and more than 10% of unaligned regions. We suggest the appropriate use of available bioinformatics resources after genome assembly to arrive at the proper bacterial identification, classification, and nomenclature to avoid erroneous species assignments and disparity due to diversity at the intraspecies level.

Characterizing the structure of styrene-maleic acid copolymer-lipid nanoparticles (SMALPs) using RAFT polymerization for membrane protein spectroscopic studies.

Membrane proteins play an important role in maintaining the structure and physiology of an organism. Despite their significance, spectroscopic studies involving membrane proteins remain challenging due to the difficulties in mimicking their native lipid bilayer environment. Membrane mimetic systems such as detergent micelles, liposomes, bicelles, nanodiscs, lipodisqs have improved the solubility and folding properties of the membrane proteins for structural studies, however, each mimetic system suffers from its own limitations. In this study, using three different lipid environments, vesicles were titrated with styrene-maleic acid (StMA) copolymer leading to a homogeneous SMALP system (∼10 nm) at a weight ratio of 1:1.5 (vesicle: StMA solution). A combination of Dynamic Light Scattering (DLS) and Transmission Electron Microscopy (TEM) was used to characterize these SMALPs. We used a controlled synthesis mechanism to synthesize StMA based block copolymers called reversible addition-fragmentation chain transfer polymerization (RAFT) SMALPs. Incorporation of the Voltage Sensor Domain of KCNQ1 (Q1-VSD) into RAFT SMALPs indicates that this is a promising application of this system to study membrane proteins using different biophysical techniques. V165C in Q1-VSD corresponding to the hydrophobic region was incorporated into the SMALP system. Continuous Wave-Electron Paramagnetic Resonance (CW-EPR) line shape analysis showed line shape broadening, exposing a lower rigid component and a faster component of the spin label.

Elucidation of stable intermediates in urea-induced unfolding pathway of human carbonic anhydrase IX.

Human carbonic anhydrase IX (CAIX) has evolved as a promising biomarker for cancer prognosis, due to its overexpression in various cancers and restricted expression in normal tissue. However, limited information is available on its biophysical behavior. The unfolding of CAIX in aqueous urea solution was studied using all-atom molecular dynamics simulation approach. The results of this study revealed a stable intermediate state along the unfolding pathway of CAIX. At intermediate concentrations of urea (2.0-4.0 M), the protein displays a native-like structure with a large population of its secondary structure and hydrophobic contacts remaining intact in addition to small confined overall motions. Beyond 4.0 M urea, the unfolding is more gradual and at 8.0 M urea the structure is largely collapsed due to the solvent effect. The hydrophobic contact analysis suggests that the contact in terminal α-helices is separated initially which propagates in the loss of contacts from centrally located ß-sheets. The reduction of 60-65% tertiary contacts in 7.0-8.0 M urea suggested the presence of residual structure in unfolded state and is confirmed with structural snap shot. Free energy landscape analysis suggested that unfolding of CAIX exists through the different intermediate states.